NASA CR- 13531 8 



CHARGE-EXCHANGE PLASMA 
GENERATED BY AN ION THRUSTER 


PREPARED FOR 
LEWIS RESEARCH CENTER 

NATIONAL AERONAUTICS AND SPACE ADMl/l^STRATION 
GRANT NSG 3038 


Annual Report 
December 1977 
Harold R. Kaufman 

Department of Mechanical Engineering 
Colorado State University 
Fort Collins, Colorado 



TECH LIBRARY KAFB, NM 


1. Report No. 

NASA CR-1 35318 

4. Title and Subtitle 


2. Government Accession No. 


CHARGE-EXCHANGE PLASMA GENERATED BY AN ION THRUSTER (U) 


7. Author(s) 


Harold R. Kaufman 


9. Performing Organization Name and Address 

Department of Mechanical Engineering 
Colorado State University 
Fort Collins, Colorado 80523 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


15. Supplementary Notes 

Grant Manager, Stanley Domitz 
NASA Lewis Research Center 
Cleveland, OH 44135 

16. Abstract 


TECH LIBRARY KAFB, NM 



3. Recipient's Catalog No. 


5. Report Date 

December 1977 

6. Performing Organization Code 


8, Performing Organization Report No. 


10. Work Unit No. 


11. Contract or Grant No. 
NSG 3038 


13. Type of Report and Period Covered 
Contractor Report 

Nov, 1, 1976 - Dec. 31, 1977 

14. Sponsoring Agency Code 


The charge-exchange plasma generated by an ion thruster has been investigated experi- 
mentally using both 5-cm and 15-cm thrusters. The results of this investigation are shown 
for wide ranges of radial distance from the thruster and angle from the beam direction. 
Considerations of test environment, as well as distance from the thruster, indicate that 
a valid simulation of a thruster on a spacecraft was obtained. A calculation procedure 
and a sample calculation of charge-exchange plasma density and saturation electron current 
density are included. 


17. Key Words (Suggested by Author(s)) 

Electrostatic thruster. Plasma physics, 
Electric propulsion 


18. Distribution Statement 

Unclassified - Unlimited 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 




’ For sale by the National Technical Information Service, Springfield. Virginia 22161 



TABLE OF CONTENTS 


Page 

INTRODUCTION 1 

APPARATUS AND PROCEDURE 3 

RESULTS AND DISCUSSION 11 

Simulation of Space Environment 11 

Survey of Charge-Exchange Plasma 12 

Radial Variation of Electron Density 15 

Angular Variation of Electron Density 17 

Saturation Electron Current Density 19 

Planar Probe Measurements 22 

Pin-Hole Probe Measurements 24 

CALCULATION PROCEDURE 27 

Sample Calculation 31 

CONCLUDING REMARKS 35 

REFERENCES 37 



INTRODUCTION 


A charge-exchange process occurs between energetic beam ions from 
an ion thruster and the escaping neutrals. As the result of this pro- 
cess, slow ions are generated. These slow ions, together with elec- 
trons that escape from the ion beam, form the charge-exchange plasma 
that surrounds an ion thruster. The presence of this charge-exchange 
plasma explains the observation of large electron currents to positive 

surfaces near thrusters.^ Without the presence of ions to neutralize 
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the electron space charge, these currents would be much smaller. The 
characteristics of this charge-exchange plasma are important for space 
missions where ion thrusters are used with high voltage solar arrays. 

The characteristics may also be of interest in predicting the environ- 
mental effects on scientific instruments that are used on an electrically 
propelled spacecraft. 

The first support period of this Grant was directed at a prelimi- 

3 

nary study of the charge-exchange plasma. A model was developed that 
gave reasonable agreement with the experimental measurements obtained 
in this study. The next support period of this Grant was directed 
primarily at the collection or direction control of the charge-exchange 
ions. Several cones, both screened and solid, were placed on the 
thruster around the ion beam. Even with the use of bias voltages, 
these cones did little more than move the effective source of charge- 
exchange ions in the downstream direction. Some estimates of electron 
collection were also made for a positive high voltage solar array on a 
spacecraft. Using an unprotected array for the beam power of an ion 
thruster, or thrusters, the collected electron current would be less 
than 10 percent of the array current below about 300 volts. Coverings 
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of insulation might permit substantially higher voltages. 

The present and final support period of this Grant was directed 
at a wider range of experimental parameters. This wider range should 
give better experimental justification for spacecraft calculations. 

The results of all support periods are summarized in this final 
report with regard to charge-exchange plasma propagation. For example, 
equations for electron density and electron saturation current density 
were obtained from the isotropic propagation model. These equations 
are the basis for scaling plasma propagation, and are therefore re- 
peated. The derivationsof these equations, however, are not repeated. 
To further aid in calculations of charge-exchange plasma propagation, 
a procedure is detailed in a separate section. 



APPARATUS AND PROCEDURE 


The vacuum facility used was the 1.2-m diameter, 4.6-m long 
chamber at the Engineering Research Center of Colorado State University. 
An 0,8-m diffusion pump maintained pressures in the mid- 10“ Torr 
range during thruster operation. A liquid-nitrogen cooled liner helped 
maintain this pressure by condensing the mercury propellent. 

Some data from earlier support periods of this Grant are included 
in this report. These earlier data were obtained with a modified 
SERT-II thruster. The modification consisted of dished molybdenum 
grids that permitted higher beam currents than the original flat grids. 

Data during the last support period of this Grant were obtained 
with a 5-cm multipole thruster designed and built for this investigation. 
This multipole thruster (see Fig. 1) had a discharge chamber 7.6 cm in 
diameter and a surrounding screened enclosure that was 17.8 cm in di- 
ameter. Tungsten wire, 0.25 mm in diameter, was used for both the 
main and neutralizer cathodes. 

The Langmuir probes used in the charge-exchange plasma survey were 
of a guarded construction, as indicated in Fig. 2. Surveys with this 
probe design extended from 26 cm upstream of the SERT-II accelerator 
grid to 49 cm downstream. Some of these earlier data are included 
herein. For the last Grant support period with the 5-cm thruster, the 
survey planes (I-VI) are indicated in Fig. 3. Langmuir probe measure- 
ments were made at several radii (from the thruster axis) in each of 
these planes. For survey planes I-III, measurements were made at 
radii of 12.3, 17.3, 18.7, 22.2, 24.0, 27.1, 29.6, 32.1, 35.4, and 
41.3 cm. For planes IV-VI the radii were 10.9, 15.2, 15.5, 19.3, 

20.7, 23.6, 26.1, 27.9, 31.6, and 37.1 cm. This survey included more 
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Stainless Steel 
Permanent Magnet 
Soft Iron 


(a) Sketch of cross section. Note screen is 
masked to 5-cm diameter. 


Fig. ^ Multipole thruster used in investigation 
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(b) Photograph of completed thruster. 
Accelerator system removed to show 
pole piece and anode construction. 

Fig. 1 Concluded. 


I 



Support End 



Fig. 2 Guarded Langmuir probe used in survey of charge-exchange plasma. 




Fig. 3 Locations of survey planes. 
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locations than the SERT-II survey, as well as extending farther up- 
stream of the thruster. A thick sheath procedure was required to 
reduce the Langmuir probe data because the probe diameter was generally 

much smaller than the Debye length. The procedure used is described by 

5 6 

Isaacson and is based on the theory and methods of Chen. 

Two additional probes were used during the 6-cm tests. These were 
the guarded planar probe (see Fig. 4.) and the "pin-hole" probe (see 
Fig. 5.). The pin-hole probe, as indicated in Fig. 5., used both a 
plain polyimide surface (Kapton, about 0.13 mm thick) and grounded 
conducting meshes attached to that surface. Both meshes were constructed 
of about 1 mm wide copper strips. The conductor behind the pin hole 
was stainless steel. Although this conductor was in close proximity 
to the polyimide, it was not cemented or glued to the polyimide. 

An ion beam target was placed at the middle of the vacuum chamber 
and electrically isolated from that chamber. All thruster operation 
was at +1000 and -500V. The 5-cm thruster was operated with a 70 
mA-equi valent neutral flow and a 50 mA ion beam. 



20 cm 




d) Probe with coarse mesh (~4cm). 



Fig. 5 Pin-hole probe configurations tested 


RESULTS AND DISCUSSION 


Simulation of Space Environment 

The vacuum chamber was negative relative to the charge-exchange 
plasma to most closely approximate the space environment. In space 
the charge-exchange ions would continue indefinitely outwards from the 
thruster. The accompanying electrons would have a drift velocity about 
equal to the ion velocity, with a much larger random velocity super- 
imposed on this drift velocity. A boundary potential that would 
collect ions and reflect electrons is thus the best simulation of 
these conditions in space. 

The negative vacuum chamber bias was established by operating the 
ion-beam target at +60V relative to the chamber. The shielding screen 
(usually grounded) surrounding the 5 cm thruster was operated at this 
same potential. Although operation with a positive bias on the target 
and shielding screen was believed to be the best simulation of space, 
the results were not found to be critical with respect to this bias. 

For example, with the target and shielding screen grounded to the 
vacuum chamber, all electron densities were found to be within a 
factor of two of the densities obtained with a positive bias. The 
reason for the lack of any large difference is believed to be the 
tendency of a plasma to become more positive than the most positive 
electrode in contact with it. Thus the reflection of electrons and 
collection of ions would normally be established. Only in very low 
density regions, where the saturation electron current is small, would 
a reflection of ions be possible. 
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Another factor affecting the validity of simulation is the back- 
ground pressure. Figure 6 shows the effect of changing vacuum chamber 
pressure by adding additional mercury flow outside the thruster. The 
lowest pressure shown in Fig. 6. was representative of normal facility 
operation, while the higher pressures required the additional mercury 
flow. The electron densities shown in Fig. 6. were from a probe 8.2 cm 
downstream of the accelerator system and 28 cm radially outwards. Ex- 
trapolation to zero pressure gives an indication of the density that 
might be expected in space. As shown by Fig. 6., this extrapolated 
density is still about 2/3 of the density at the normal operating pres- 
sure. The accuracy of this extrapolation is limited by the estimated 
+ 50 percent accuracy of the probe density measurements. On the other 
hand, some of background pressure during normal operation must be due 
to gases other than mercury. The charge-exchange cross sections of 
mercury beam ions with these other gases should be much smaller than 
with mercury. A decrease in pressure due to removal of these other 
gases should therefore not have a significant effect on the charge- 
exchange plasma density. In short, most of the charge-exchange plasma 
around the thruster would still be expected to be present in space. 

Survey of Charge-Exchange Plasma 

The charge-exchange plasma surrounding the 5-cm ion thruster was 
surveyed using Langmuir probes. The electron density was the most 
reproducible plasma property, as well as being the one of most interest. 
A survey of this density is shown in Fig. 7. This survey is similar to 
those obtained earlier for the 15-cm thruster, except that a smaller 
ion-beam diameter was used and the survey extended further upstream 
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of the thruster. Both of these changes permitted a more rigorous 
evaluation of the transport model and its validity for spacecraft 
calculations. 

The electron temperature was also obtained from the Langmuir 
probes and ranged from about 3 to 4 eV in the charge-exchange region. 
This was slightly higher than the 2.5 eV average obtained for the 
same region with the 15 cm thruster. A refractory metal neutralizer 
was used with the 5-cm thruster, while a plasma bridge neutralizer 
was used with the 15-cm thruster. The electron temperature difference 
was probably due to the size or neutralizer difference for the two 
thrusters. 


Radial Variation of Electron Density 


The radial variation of electron density is, according to the 
isotropic transport model, an inverse-square relationship. In this 
model the production of charge-exchange ions is calculated by approxi- 
mating the neutral density downstream of the thruster with the known 
solution downstream of a circular sharp-edged orifice. For the purpose 
of simplifying the calculations and obtaining a conservative answer, 
the ion beam was assumed to be concentrated near the axis of the 
neutral efflux. With the further assumptions of plasma neutrality 

g 

and that the charge-exchange ions leave at Bohm velocity, the elec- 
tron density and saturation electron current density are 


1.49 X 10^^ (1-n,) o,, A 

>-b \ -Vf; 


S 


( 1 ) 
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J6 


^ 2.62 X 10- (l-n,) A ^ 


^ 


( 2 ) 


_3 

where electron density n„ is in m , electron current density j is 

el '"e 

2 

in A/m^ beam current J|^ is in A, charge-exchange cross section 
2 

is in m^ atomic weight A is in amu, beam radius rj^ is in m, radius 

from charge-exchange ion source R is in m, neutral and beam electron 

temperatures and are in °K, and propellant utilization is 

-19 2 

dimensionless. Substitution of 200.6 for A, 6 x 10 m for 
9-11 

°ce’ ^ typical value of 500°K for T^, and a typical observation 
(in this investigation) of 58,000°K (5eV) for T^, the preceding 
equations become 

ng = 3.3 X 10^^ 


j- = 0.14 J, 




(4) 
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It should be noted that lower values of than 5eV have been observed, 
but the value used here is consistent with the usual small thruster injec- 
tion voltage of about 20 and the typical ratio of 0.3 for electron temperature 

13-14 

divided by injection voltage. It should also be noted that elec- 

tron temperature does not appear in Eq. (2), hence the exact value of 
this variable should not significantly affect observed values of sat- 
uration electron current. 

From previous work under this Grant, this isotropic model should 
be compared to experimental plasma properties normal to the ion beam 
direction. Rearranging Eq. (4), we obtain 


"e '"b \ 


/ (l-n ) = 3.3 
o u 


X 10^^ 


(5) 
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A curve for this equation is shown in Fig. 8., together with experi- 

mental values of the parameter n^ rj^ n^j / The data points 

shown in Fig. 8., are based upon individual measurements of electron 

density (not faired values) and range from 86 to 94 degrees relative 

2 

to the ion beam direction. The data points clearly follow the 1/R 

trend of Eg. (5). Further, the 5-cm data obtained during this last 

Grant support period cover a much wider range of R/r^^ than the 

preceding 15-cm data. From Fig. 8., then, it appears justified to 
2 

use the 1/R relationship for spacecraft calculations. 

A final point should be made about Fig. 8. concerning the assumed 
point of origin for the charge-exchange plasma. From the data obtained 
with the 15-cm thruster, the origin was assumed previously to be one 
ion beam radius downstream of the accelerator system. A more recent 
evaluation, which includes the 5-dm data, indicates that one radius of 
the screen surrounding the thruster would be a better choice than one 
ion-beam radius. This origin change is only about 4.5 cm for the 15-cm 
thruster, but results in considerably better agreement between 15-cm 
and 5-cm data. This redefinition of the origin point is also con- 
sistent with results observed using cones downstream of the ion 

4 

thruster. 


Angular Variation of Electron Density 

Having justified the inverse-square relationship of electron 
density with radial distance, it is appropriate to examine the angular 
variation. The appropriate electron-density parameter can be obtained 
by rearranging Eq. (3). 



Radius ratio, R / 


Fig. 8 Effect of radius ratio R/rj^ on electron density. The 
electron density is in a corrected or normalized form 
to permit inclusion of data from both the 5-cm and 
15-cm thrusters. The angle from the beam direction is 
from 86 to 94 degrees for the data shown. 
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The electron density parameter n^ rj^ R plotted 

against the angle relative to beam direction in Fig. 9. As in Fig. 8, 
actual electron density measurements were used, not faired values. 

The point of origin for the charge exchange ions was again assumed to 
be downstream of the accelerator system by one radius of the screen 
surrounding the thruster. The average level of the data near 90 
degrees is consistent with Eq. (6), as would be expected from Fig. 8. 

The data for angles less than 30 degrees were omitted because they 
approached, or were in, the ion-beam volume. Between 30 and 90 degrees 
the data spread is substantial. This spread was observed earlier and 
is believed due to the particular beam profiles and angular divergence 

3 

obtained from the thrusters and operating conditions. The data 
spread at angles greater than 90 degrees is about equal to that at 90 
degrees, indicating that the radial variation at angles greater than 
90 degrees would be similar to that shown in Fig. 8. The 15 cm data 
tend to fall below the 5 cm data at angles greater than 140 degrees, 
but the total data spread is still believed to approximate the 
accuracy of the electron density measurements. In view of the larger 
ranges of experimental variables covered in Fig. 9, the experimental 
variation with angle shown in Fig. 9 is preferred over the semiemperical 
equation developed earlier. 

Saturation Electron Current Density 

For the isotropic model, the saturation electron current density 
is given by Eq. (4). Rearranging this equation to obtain the electron 


Electron density parameter, 
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30 60 90 120 ISO 180 

Angle from beam direction, degrees 


Fig. 9 Variation of electron density parameter with 
angle from beam direction. Assumed point of 
origin is one thruster (ground screen) radius 
downstream of the accelerator grid. 
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current parameter, we obtain 

Je '■b \ ' '^b <’-%> = '>•’'* ’ (?) 

where is the current density (A/m ) for any surface that is 
sufficiently positive to collect all the arriving electrons. The 
most extensive data were obtained in terms of electron density rather 
than saturation electron current density. To convert from electron 
density to current density, we use 

J'e " "e ^e ^ ^ 

1 

where v is the average electron velocity, (SkT^/nm^)^. With the 

substitution of numerical values, Eq. (8) becomes 

jg = 2.49 X lO"''^ Tg^ . (9) 

The particular temperature used in this equation needs a word of 
explanation. The electron temperature in the ion beam has been found 
to be about twice the electron temperature in the charge exchange 
region. The 58,000°K (5eV) temperature used to obtain Eq. (3) was 
measured in the ion beam. The corresponding electron temperature in 
the surrounding charge-exchange plasma was about 29,000°K (2.5eV). 
Inasmuch as we are concerned with the charge exchange region, the 
temperature to use in Eq. (9) should be for that region. That is, 
about 2.5eV for the 15 cm data and 3.5eV for the 5 cm data. Note that 
the electron temperature is required because we started with electron 
density measurements. As shown by Eq. (2), the electron temperature 
is not a variable in the collection of electron current if the current 
is measured directly. 
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Equation (9) was used to convert the electron densities to 
equivalent electron current densities. The corresponding electron 
current parameters were plotted in Fig. 10 against the angle from 
the beam direction. The data distribution in Fig. 10 is similar to 
that shown for electron density in Fig. 9. At angles less than 90 
degrees from the beam direction the data spread is large due to the 
influence of the particular beam profile. At angles of 90 degrees, 
or more, the data spread approximates the accuracy uncertainty of 
Langmuir probe measurements. 

Planar Probe Measurements 

The electron current densities of Fig. 10 were obtained from 
electron densities and electron temperatures, both of which were 
obtained from Langmuir probes. Because the probe diameter was small 
compared to the Debye length, there was no direct measurement of the 
saturation electron density. The planar probe was designed to make 
such a direct measurement. The Debye length in the vicinity of the 
planar probe was 2-3 cm. The 5 cm guard around the probe was therefore 
sufficient to eliminate edge effects around the central 10 cm square 
probe surface (see Fig. 4). The potential of the planar probe was 15 
to 25 V above the local plasma potential, which was determined from 
Langmuir probe survey. The potential should therefore have been suf- 
ficiently positive to assure collection of all arriving electrons. 

The current collected by the planar probe was 0.22 mA. For the 
radial location used, together with thruster operating conditions, 
this current corresponded to an electron current parameter of 0.050. 



Electron current parameter, 
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Angle from beam direction , degrees 

Fig. 10 Variation of electric current (density) 
parameter with angle from beam direction 
Assumed point of origin is one thruster 
(ground screen) radius downstream of 
accelerator grid. 
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This value, at an angle of 122 degrees from the beam direction, is 
near the center of the data spread in Fig. 10. The planar probe 
measurement thus constitutes a separate and independent check on the 
calculated current densities used in Fig. 10. 

Pin-Hole Probe Measurements 

The major objective of the overall investigation was to evaluate 
the charge-exchange plasma surrounding an electric thruster. As such, 
the pin-hole probe was considered a minor experiment that could be 
added with only a small increase in the total experimental effort, 
inasmuch as the charge-exchange plasma was already available. 

The electron density in the vicinity of the plasma probe was about 
10 -3 

9 X 10 m . This density, together with the 3.5eV electron tempera- 
ture, resulted in a Debye length of about 5 cm. Comparing this length 
to the probe geometries investigated, the plain probe (Fig. 4 (a)) was 
larger than the Debye length, the coarse mesh (Fig. 4 (b)) was roughly 
equal to the Debye length, and the fine mesh was smaller than the 
Debye length. The purpose of the mesh was, of course, to interrupt 
either surface conduction or plasma-sheath focusing, thereby reducing 
the current collected by the pin hole. 

The current-collection data for the pin-hole probe are shown in 
Fig. 11. The current at a given collection voltage shows a systematic 
decrease from the plain probe to the coarse mesh to the fine mesh. The 
magnitude of the decrease, though, is much smaller than the area re- 
duction from the total probe to the area enclosed by one mesh square. A 
shortcoming of the data shown in Fig. 11 is that collection voltages 
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Fig. n Current conected by pin-hole probe. 

(See Fig. 5 for probe configruations. ) 




-26- 


were not high enough to encounter the sharp increase in collection 
current that is normally found. The effect of insulator area is 
usually larger above this increase, so that the three configurations 
of Fig. 11 might be expected to diverge more at higher voltages. 


I 



CALCULATION PROCEDURE 


The assumed point of origin for the charge-exchange plasma surround- 
ing a thruster is indicated in Fig. 12. For a single thruster this point 
is one thruster radius r^ downstream of the accelerator grid and on the 
thruster centerline, as shown in Fig. 12(a). For a circular array of 
thrusters (Fig. 12(b)), the effective point of origin would be expected 
(from conical collector tests^) to be one radius of the thruster array r| 
downstream of the array center. For a noncircular array, the effective 
array radius r| would be expected to be between the maximum and minimum 
values from the center of the array. A noncircular array would probably 
also cause some departures from axial symmetry. But the magnitude of 
these departures cannot be estimated at the present time, except that the 
charge-exchange plasma would probably depart less from axial symmetry than 
the array. 

The charge-exchange plasma at some point of interest (see Fig. 12(a)) 

is a function of the radial distance from the assumed point of origin, R, 

and the angle relative to the beam direction, e, in addition to the 

thruster parameters of beam radius rj^, beam current Jj^, and thruster 

propellant utilization As shown by Figs. 9 and 10, the electron 

★ 

density n^ and saturation electron current density are effectively 

correlated by parameters including these variables. Rather than use 


* 


This current density is the value that would be obtained at any surface 
that is sufficiently positive to collect all electrons approaching it 
and sufficiently large to make edge effects negligible. 
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^ — Point of interest 


R 



(a) Single thruster 




I 


I 
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(b) Thruster array 


^Assumed point 
of Origin 


Fig. 12 Parameters for calculation of 

charge-exchange plasma around thruster. 
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Figures 9 and 10 for calculations, the upper envelopes of the data dis- 
tributions in these figures are tabulated in Tables I and II as functions 
of the angle e. To calculate the electron density and saturation electron 
current density, select values from Tables I and II for the angle 8 of 
interest and substitute in the equations below: 

2 

n^ = [J|^ ^ *^b *^^\^tTable I Param] m"^ (10) 

.I'e = 0-n^,) / rj^R\][Table II Param] A/m^ (11) 

The beam current is in A, the propellant utilization is dimensionless, 
and both the beam radius r. and the radius R are in m. 
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Table I - Electron density parameter as a function of 
angle from beam direction. (From Fig. 9.) 


e, deg. 

Param. 

e, deg. 

Param. 

e, deg. 

Param. 

30 

2.3x10^^ 

80 

6.6x10^^ 

130 

12 

l.Oxlo'^ 

40 

1.8x10^^ 

90 

12 

5.0x10“^ 

140 

6.2x10^^ 

50 

13 

1.4x10*'^ 

100 

12 

3.5xl0'^ 

150 

3.7x10^^ 

60 

13 

1.1x10 

no 

12 

2.4x10 

160 

2.2x10^^ 

70 

12 

8.6xl0'^ 

120 

1.6x10^^ 

170 

1.3x10^^ 

Table II - Electron current density parameter as a 
function of angle from beam direction. (From Fig. 10 

e, deg. 

Param. 

e, deg. 

Param. 

6, deg. 

Param. 

30 

1.2 

80 

0.32 

130 

0.052 

40 

0.94 

90 

0.23 

140 

0.032 

50 

0.72 

100 

0.17 

150 

0.019 

60 

0.55 

no 

0.12 

160 

0.011 

70 

0.42 

120 

0.082 

170 

0.064 



Sample Calculation 


The electron density and saturation electron current density have 
been calculated for an 8-cm thruster as a sample calculation. The 
operating conditions selected were obtained from literature^^’ and 
are: 

Total neutral flow, 0.0914 A-equi valent 
Net accelerating potential difference, 1205. V 
Beam current, 0.072 A 
Propellant utilization, n^j, 0.79 
The pertinent physical dimensions are: 

Radius ion beam, rj^, 0.04 m 
Radius of ground screen, 0.084 m 

The operating conditions given assume the use of a small hole accelerator 
grid, which permits higher propellant utilization than would otherwise be 
practical . 

The point of interest for these calculations was arbitrarily assumed 
to be 1.0 m radially out from the thruster centerline and 0.5 m upstream 
of the accelerator grid. The assumed point of origin for the charge-exchange 
plasma is one ground-screen radius downstream of the accelerator grid (see 
Fig. 12(a)). This gives a radius R of 

R = (1.0^ + 0.584^)^^"^= 1.16 m. 

For the point of interest selected, the angle relative to the beam direction 
is 

0 = 90 + arc tan (0.584 / 1.0) 


= 120 °. 
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1 9 

From Table I, the electron density parameter for 120 degrees Is 1.6 x 10 . 

(Interpolation could have been used for an angle between those given in 
the table.) Substitution of the numerical values for the various parameters 
in Eq.(lO) yields an electron density of 

n = [0.072^(1-0.79) / 0.04 x 1.16^ x 0.79][1.6 x 10^^] 

10 -3 

= 4.1 X lO'^ m 

From Table II, the electron current density parameter for 120 degrees is 
0.083. Substitution in Eq. (11) then yields a saturation electron current 
density of 

j = [0.072^(1-0.79) / 0.04 x 1.16^ x 0.79][0.082] 
e 

= 0.0021 A/m^. 

This current density is, of course, the value that would be collected by 
any surface that was sufficiently positive to collect all electrons directed 
toward that surface. It is also assumed that the surface in question is 
approximately normal to the radius R. 

The calculated values for n^ and j’ represent the upper edge of the 
distributions shown in Figs. 9 and 10, hence should be conservative (high) 
values. From the distributions shown near the angle of 121 degrees, the 
actual values may be as much as a factor of 3 lower. 

The foregoing is true for electron temperatures close to those 
observed in the experimental investigation. For the current density j^, 
the result should be unaffected by changes in electron temperature. This 
is because Eq. (11) is based on Eq. (2), which has no electron temperature 
in it. 
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For the electron density n^, an effect of electron temperature would 
be expected. As shown by Eq. (1), upon which Eq.(lO) is based, the elec- 
tron density should vary inversely as the square root of electron tempera- 
ture. The measurements of this investigation indicated 2.5 and 3.5 eV for 
the 15-cm and 5-cm thrusters, respectively. The electron temperature 
within a 30-cm thruster beam was found to be 0.35 eV. Including the 
observation in this investigation that beam electron temperature is about 
twice the electron temperature in the surrounding charge-exchange region, 
there is an uncertainty in electron temperature of up to a factor of 20. 
This temperature uncertainty translates into a factor of 4 or 5 in elec- 
tron density uncertainty. Specifically, n^ could be 4 or 5 times as high 
as the calculated value due to a lower electron temperature. 

Although the total electron density uncertainty due to electron 

temperature must be of the order of this factor of 4 or 5, the probable 

uncertainty in electron density is believed to be smaller. The observed 

electron temperature appears to decrease with increasing thruster size. 

Such a trend has also been observed elsewhere.^^ A similiar decrease in 

Maxwellian electron temperature has also been observed within the dis- 

1 8 

charge chamber for increasing thruster size. Inasmuch as the 8-cm 
thruster size in the sample problem is between the 5-cm and 15-cm sizes 
used in this investigation, the electron temperature might be expected 
to also be intermediate. Conversely, if a calculation were being made 
for a 30-cm thruster, a more serious consideration should probably be 
given to the effect of a lower electron temperature. 


- 34 - 


The sample calculation presented above is for a single isolated 
thruster. An array of thrusters will produce a total charge-exchange 
plasma that is greater than the sum of that produced by the same number 
of isolated thrusters. This is because the beam of each thruster passes 
through the neutral efflux of all the thrusters, thus increasing the 
charge-exchange ion production over the value that would be obtained 
from an isolated thruster. Including this interaction requires tedious, 
but straightforward free-molecular- flow calculations. Some simplifica- 
tion of these calculations is possible by asssuming that the beam from 
each thruster is concentrated near the axis of that thruster. 



CONCLUDING REMARKS 


The validity of simulating the charge-exchange plasma generated by a 
thruster in space was investigated by varying the mercury background 
pressure in the vacuum chamber. The resulting measurements indicate that 
most of the charge-exchange plasma found in the vacuum chamber would also 
be present in space. 

The survey of the charge-exchange plasma surrounding a 5-cm thruster 
gave results generally in agreement with earlier surveys using a 15-cm 
thruster. 

The variation of electron density was investigated over a wide range 
of radial distance relative to the beam radius. The data for this ex- 
tended range agreed closely with the inverse-square theoretical model, and 
should justify the calculations at even larger radial distances for space- 
craft environments. 

The angular variation of electron density and saturation electron 
current density in the charge-exchange plasma was found to follow a slightly 
different trend than the semiemperical equation obtained earlier. Plots 
of the experimental data appear to be the best basis for predicting angular 
variations of these plasma parameters. These plots include data sets from 
both the 5-cm and 15-cm thrusters, with excellent agreement between the 
two sets. 

Electron temperature is a variable that will affect the electron 
density, but not the saturation electron current density. Because the 
effect is through the square-root temperature, the estimates of electron 
density should still be useful. 
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The experimental results are also presented in tabular form to aid 
calculation. A calculation procedure is included, together with a sample 
calculation of electron density and saturation electron current density. 
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